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We formulate a QCD sum rule to find the three momentum dependence of the peak position of 
the vector meson spectral density in nuclear medium. We find less than 2 % (0.1 %) shift of the 
peak position at nuclear matter density and at (q = 0.5GeV/c) for the p, uj (<j>) mesons. However, 
at higher density and three momentum, its effect becomes non- negligible and the relevance of our 
result in the dilepton spectrum in A-A and p-A reaction is discussed. 



The properties of vector mesons in nuclear medium have attracted a lot of interest because of their potential role 
to experimentally observe a non-perturbative aspect of QCD, namely the restoration of spontaneously broken chiral 
, symmetry at finite temperature or density, through dileptons from A-A or p-A reaction jij . 

Many model calculations have been performed to calculate the vector meson mass shift at finite density By now, 
there seems to be a consensus that the average peak position of the vector meson spectral density at zero three 
' momentum (q = 0) will shift down at finite density . 

Indeed, there already exist dilepton data from the CERES collaboration, which report an enhancement of low 
mass dileptons below the p meson invariant mass in S-S and recently from Pb-Pb collisions at CERN Q. So far, all 
conventional collision models failed to explain the enhancement except when the p meson mass is allowed to decrease 
in medium as predicted by theoretical calculations ||. However, before coming to any definite conclusions, it is 
necessary to consider all possible conventional mechanisms. 

In nuclear medium, in addition to possible change in the vector meson mass, there will be breaking of Lorentz 
t^J- , invariance and hence two independent polarization directions of the vector mesons. Each polarization will have a 
different dispersion relation, which to leading order in the three momentum would be modified to ui 2 — (l+a)q 2 — my = 
with a ^ 0. Suppose, experimentally, one detects a dilepton with energy u> and three momentum q, then the average 
peak due to the vector meson will appear at M 2 = (m v +aq 2 ), so that if a < 0, the strength will be shifted downwards 
£^ , for dilepton pairs with q ^ 0. Hence, it is important to estimate the finite q effect. 

In this letter, we formulate a QCD sum rule to find the finite q effect to leading order in density for the transverse 
and longitudinal direction of the p,u> and <fi mesons. This formalism also provides the first attempt to estimate the 
^ ' leading q dependence of the V — N T matrix, when it has a small off-shell dependence. 

Let us consider the correlation function between the vector current with p, w and (f> meson quantum numbers, 
JP' U — ^(m7 m u =f dj^d) and = sj^s in the nuclear medium. 
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n^w.q) - * / d A xe* x {T[J^x)J v {Q)]) n . m .. (1) 

Here () n .m. denotes the nuclear matter expectation value. In general, because the vector current is conserved, the 
correlation function in eq.(^) will have two invariant functions H. 

n MV (u, q) - n T q 2 Pl„ + n L q 2 P L ^, (2) 

where for q = (u),q) and medium at rest, we have, Pq = Pq, = P^ = 0, Pfj = Sij — q^q^/q 2 , and P^ v — 
(q^q u /q 2 — g^u — Pp^)- In the limit when q —> 0, there is only one invariant function Hi(w,0) = Ht{^,0). In this 
work, we will formulate a QCD sum rules for both IIf,(a;,q) and Hr(o;,q) at finite q. 

The starting point is the energy dispersion relation at finite q. For small q 2 < ui 2 , we can make a Taylor expansion 
of the correlation function such that, 

Ren L , T (w 2 , q 2 ) = Re (U° L jT (w 2 , 0) + n*^ 2 , 0) q 2 + ■•) 
' 2 ( p(u,0)l T p(u,0)l T \ 
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here p(u, q) = l/nlmIl R (u, q), and R denotes the retarded correlation function. As we will discuss below, the left 
hand side (Ren) is known only to leading order in q 2 and since we are interested in the leading behavior anyways, 
we will look at the dispersion relation for n 1 in eq.(|^). 
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The real part of eq. (gj) is calculated via the Operator Product Expansion (OPE) at large 
The full polarization tensor will have the following form. 



oo with finite q. 



rV(^,q)= (q^-g^) -coln|Q 2 | + £ d c -^A d > d {i 



+ £ S , T=2 [ c d,r9^ 11 ■ -q^A^^Jn.m.) + •■] , (4) 

where, Q 2 — q 2 — to 2 . Here, A d,T (n.m.) represents the nuclear matter expectation value of an operator of dimension 
d and twist r = d — s, where s is the number of spin index. These operators are defined at the scale Q 2 and the c's 
are the dimensionless Wilson coefficients with the running coupling constant. This way of including the density effect 
is consistent at low energy pof . The first set of terms in eq.(|^) come from the OPE of scalar operators, the second 
set from operators with non zero spin s. 

To linear order in density, the matrix elements are related to the nucleon expectation values of the operator via 

A(n.m.) = A Q + n n A N , (5) 

where Aq is the corresponding vacuum expectation value, A n the nucleon expectation value and n n the nuclear 
density. 

As in the vacuum, we will truncate our OPE up to dimension 6 operators. This implies that in our OPE in eq.(Q), 
we will have contributions from (r, s) — (2, 2), (2, 4), (4, 2). The nucleon matrix elements of the t = 2 operators are 
well known. The r = 4 matrix element appearing in the p, lo sum rule are similar to those appearing in electron D1S 
jnj and have been estimated jl2|,|f3| up to about ± 30% uncertainty from available D1S data from CERN and Slac. 

The q dependence coming from the first line of eq.(^), namely the contribution from the scalar operators, come from 
the q dependence in Q 2 . These form the so called "trivial" q dependence, and comes from replacing lo 2 — ► lo 2 — q 2 
when going from zero to finite three momentum. Here, we are not interested in these trivial dependence and also not 
in the possible change in the scalar mass my. Consequently we do not need the nucleon expectation value of the 
scalar operators. Operators with spin also partly contribute to the change in the scalar mass and hence also to the 
trivial changes. However, these spin parts also give the non-trivial q dependence. A prescription to find the nontrivial 
q 2 dependence in the OPE is to first calculate the total q 2 term n 1 and then subtract out the trivial dependence. 

+ _> (d (6) 

L0 n lo z 2 L0 n+Z 

Following this prescription, we find the following contributions from the r — 2,4 operators. 

ni, r H/p„ = ^ + (7) 

For p, lo, the transverse (T) and longitudinal (L) parts give, 



&2 = i\cl, 2 - \cl 2 )mA u 2 +d + (C G 2 - C G 2 )mA G (M 



bl = {\CI A - \cl A WA\+ d + ( 9 -C? A - 5C G 4 )m 3 A 4 G + \m (-(1 + 0){K X + *-K 2 + ^K») + ± 1) 



b2 = -\c% 2 mAl +d - Cl 2 mA 



(9) 
(10) 



= (^L - -Cl A WAT- d + (C G 4 - 5C G 4 )m 3 A 4 G + -(f + 0) (tf 2 - ) . (Ill 
where ± refers to the p and lo case. Here, m is the nucleon mass and for even n, 



A« n = 2 [ to"" 1 [<?(:£, Q 2 )+g(a;,Q 2 )], 
Jo 

A G = 2 f dxx n ~ l G{x,Q 2 ), 
Jo 



(12) 



where q(x,Q 2 ) and G(x,Q 2 ) are the quark and gluon distribution functions. We will use the HO parameterization 
for these obtained in ref [Ful which should be used with the Wilson coefficients C's in the MS scheme [n5f. Terms 
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proportional to K's come from r — 4, s — 2. We will use the set of K values obtained in ref. jl2); ( K l , K 2 , K\ d , K 9 ) = 
(-0.173GeV 2 ,0.203GeV 2 ,-0.083GeV 2 ,-0.238GeV 2 ) and we take (3 = 0.5. For the r = 4 operators, we neglect the 
Q 2 dependence. 
For the <fi meson, 

b\ = (2C| >2 - 1Cl^mA\ + (2C G 2 - 2Cl 2 )mA§ (13) 
bl = (9C|, 4 - 10C« A )m 3 A s 4 + (9C G 4 - 10C G 4 )m 3 A G - ^mC, (14) 
^ = -2Cl )2 m^ - 2C G 2 mA G (15) 
& 3 = ( 2 Cf,4 - 10C« 4 )m 3 A^ + (2C* 2 G 4 - 10C G , 4 )m 3 A G - 9mC, (16) 

where C is defined from (iV(p)|sD p D„TO s s|iV(p)) = (p M p„ - \m 2 g llv )C. We let C = ~xJN\sm s s\N) = -0.1548 • 
a; s GeV 2 . Where we used the same number for the strange content of the nucleon as in ref. |l6| with the normalization 
of (p\p) = 2u>(2ir) 3 5 3 (p — p). x s is an unknown number, however, assuming its order to be similar to the ratio j4|/A|, 
we will take x s = 0.04. 

Table 1 summarizes each contribution to the b coefficients at Q 2 — lGeV 2 . As can be seen from the table, for 
the transverse part of the p, uj 1 the well known twist-2 quark contribution dominates both the 62,^3- However, for 
the longitudinal part, both the quark and gluon twist-2 part contributes at order a s and the twist-4 part becomes 
important, which has larger uncertainty. This is also roughly true for the </> meson. At higher Q 2 values, the values 
of 6's decreases in general. 

In the vacuum, the spectral density appearing in the left hand side of eq.(^|) is modeled with a pole and a continuum: 
8ir 2 p(u) = FS(u 2 — m 2 /) + cq8(u 2 — So). In our case, we allow the three parameters to vary non-trivially by a term 
proportional to q 2 and the nuclear density n n . Such as, F — » F + f ■ q 2 , my — > rriy + a ■ q 2 , So — > So + s ■ q 2 . Also, 
when using the dispersion relation in eq.(3), there is a possible ambiguity in the subtraction constants and we have 
to know if there exist singularities in the spectral density of the form 8(u 2 ) or 8'{u 2 ). The contribution proportional 
to 5{u 2 ) is not known. However, the other singularity can be unambiguously calculated from the the nucleon-hole 
contribution. This term is called the scattering term || and gives a non-trivial contribution to the longitudinal 
direction. Hence, the spectral density of n 1 (w) is, 

8ttV(u) = / • 8{u 2 - my) - a ■ FS'(u 2 - m v ) (17) 
-s ■ c S(u 2 - So) + 87r 2 n n b scatt 5'(u 2 ), (18) 

where, b sca tt — l/(4m) for the longitudinal p, uj and zero otherwise. For conventional reason, 47r 2 will be used for <fi 
meson instead of 87r 2 . 

The Borel sum rule with the 8(u 2 ) ambiguity subtracted out is obtained by taking the Borel transform of u^Ii 1 ^). 
This gives, 

(m 2 p f + F(l - g|)a) e-<'^ - coS e-^ M ' s 

2 ( & 2 



= ^^{M2~2m~ bacatt )- (19) 

The Q 2 dependence in the 6's coming from the twist-2 operators changes to the M 2 dependence |17[| . The vac- 
uum parameters are first determined from the vacuum sum rules with parameters given as in ref. juT This gives 
(m v ,F,S ) = (0.77 2 ,1.48,1.43)GeV 2 for the p,w and (m 2 ,F,S ) = (1.02 2 , 2.19, 2.1)GeV 2 for the (j>. Then the pa- 
rameters are determined by least square fit method. The Borel interval for the transverse p meson is determined by 
requiring that the contribution from dim 6 operators are less than 35% of the dim 4 contributions, which determines 
the M^ in ~ lGeV 2 . The maximum Borel mass is determined by requiring that the continuum contribution is again 
less than 35% of the first term in the OPE, which gives M^ ax ~ 2.3GeV 2 . But it turns out that even if we choose 
higher M^ ax , the result differ by less than 10%. The uncertainty here comes from the less reliable number of the 
6|(r = 4), which gives less than 10 % uncertainty in the final answer. Overall, combining the uncertainty coming from 
the continuum, we expect 40% uncertainty for the numbers for the transverse p and also for the transverse lu, <fr, which 
can be analysed in a similar fashion. For the longitudinal directions of p,u>, the dim 4 operators are suppressed by a s 
and the dim 6 operators contribute at tree level. So here we choose the M m i n ~ 2.5GeV 2 , at which the dim 4 and 
dim 6 contributions are of similar order and take M max ~ 3.5 to 5GeV 2 . Here, we expect 50 % uncertainty. Similar 
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analysis can be done for the longitudinal (f>. The results are shown in table 2, the numbers are all at the nuclear 
matter density. For higher density, one can just multiply the numbers with the relevant ratio to the nuclear matter 
density. Our results are consistent with effective hadronic model calculations. First, as can be seen from table 2, the 
scattering term decreases the longitudinal a value. This is consistent with the known result from the nucleon-hole 
contribution fl6|. Second, recent calculation using resonance-nucleon hole contribution for the transverse part shows 
an attraction Eq]. This is also consistent with our a values for the transverse part of the p,uj. 

As discussed before, a non- vanishing a will shift the average peak position by AAI = U m v + aq 2 — my, even 
if there is no change in the scalar mass my- With the values of a obtained, we have plotted the fractional change 
AM/my for the p meson in Fig. 1 a). The solid lines denote the result at nuclear matter density, and the dashed 
lines that at 3 times nuclear matter density. The results for the u> meson look similar to Fig.l a). Fig 1 b). shows the 
result for the ip meson. As can be seen from the figures, at nuclear matter density, the expected shift is small, and will 
not wash out any shift of the peak position expected from the change of the scalar mass But, at higher densities 
and higher q ~ lGeV, its effect becomes quite large. In the dilepton production, the relative contribution from the 
transverse and the longitudinal direction depends on the angle between the total and relative three momentum of the 
outcoming dileptons |is|| . However, the transverse direction has always a larger contribution. In addition, our result 
show that the residue (F + /q 2 ) increases more for the transverse direction as q increases. Therefore, the fact that 
a < for the transverse direction will shift the dilepton spectrum downwards for finite q. In addition, the contribution 
of the r = 2 matrix element is the same for the axial vector axial vector correlation function as in the case of vector 
and vector. So the dominating effect that pushes a negative has nothing to do with chiral symmetry breaking and 
its restoration. It is not clear yet, how much this effect will be in an actual RHIC. For that, one has to do a model 
calculation and such kind of calculation is important to really pin down the possible change of vector meson mass at 
finite density. 

We would like to thank Bengt Friman for introducing me to the importance of finite q effect, which initiated this 
study. We thank T. Hatsuda and H.C. Kim for useful comments. This work was supported in part by the Korean 
Ministry of Education through grant no. BSRI-96-2425, by the KOSEF through grant no. 971-0204-017-2 and through 
the CTP at Seoul National University. 

Table 1. 62(9) (&2(G)) represents the contribution of quark (Gluon) operators to 62. 632 (^34) represents the contri- 
bution of the r = 2 (r = 4) operators to 63. The units for 62, 63 are GeV and GeV 3 respectively 
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bh 
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&J(G) 


&I2 


O34 
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-0.021 


-0.014 


0.047 


0.096 


0.417 


-0.025 


0.242 


0.060 


OJ 


-0.021 


-0.014 


0.047 


0.096 


0.417 


-0.025 


0.242 


0.138 


4> 


-0.004 


-0.028 


-0.003 


0.052 


0.074 


-0.049 


-0.002 


0.015 



Table 2. Results for the parameters at nuclear matter density. The values are from best fit of the Borel sum rule in 
eq.(|l^). The values in the bracket are the results without the scattering term. 
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Borel Interval GeV 



Transverse p 
Transverse ui 
Longitudinal p, u 

Transverse (f> 
Longitudinal <f> 



-0.065 
-0.040 
0.021 
(0.061) 
0.004 
0.009 



0.137 
0.120 
0.068 

(-0.042) 
0.010 

-0.001 



-0.008 
0.009 
0.027 

(0.042) 
0.009 
0.009 



1 <- 
1.3 
2.5 

0.9 
2.0 



2.3 

- 2.5 

- 3.5 

- 2.0 

- 3.0 
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FIG. 1. a) The fractional change AM/mv of the peak position of the p as a function of q. The solid (dashed) lines show the 
results at nuclear matter (three times nuclear matter) density. The positive changes correspond to the longitudinal direction, 
the negative changes correspond the transverse directions, b) The fractional changes of the as a function of q. The solid 
(dashed) lines show the result at nuclear matter (three times nuclear matter) density. The larger changes correspond to the 
longitudinal direction, the smaller changes correspond the transverse directions. 
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